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Quintessence theories for cosmic acceleration imbue dark energy with a non-trivial dynamics that 
offers hope in distinguishing the physical origin of this component. We review quintessence models 
with an emphasis on the dynamics and discuss classifications of the different physical behaviors. 
The pros and cons of various paramctcrizations are examined as well as the extension from scalar 
fields to other modifications of the Friedmann expansion equation. New results on the ability of 
cosmological data to distinguish among and between thawing and freezing fields are presented. 



I. INTRODUCTION 

Understanding the acceleration of the cosmic expan- 
sion is a landmark problem in physics, impacting gravita- 
tion, high energy and quantum physics, and astrophysics, 
and likely to revolutionize one or more of these fields. The 
direction in which to look for a solution is almost wholly 
unknown currently. Though there is no shortage of sug- 
gestions, most are far from a first principles explanation 
of how such physics arises. 

Perhaps the simplest propos al - Einstein's cosmologi- 
cal constant A ( Einsteir]Jl917l ) - is correct, though even 
so we have as yet no understanding of why it would 
arise, with the magnitude needed to explain accelera- 
tion occurring near the present epoch. That puzzl ement 
can be broken into two severe p roblem s (|CarrollL l200lt 
iPadmanabhanL 120031 : IWeinbereL 1 1989( 1: the fine tuning 
problem of how A appears with a magnitude (energy den- 
sity or energy scale) so far from the natural (Planck) scale 
defined by fundamental constants, and the coincidence 
problem of why acceleration appears in our recent past, 
at a cosmic scale factor within 2 of the present value out 
of perhaps 10 28 since inflation. The cosmological con- 
stant is addre s sed in far great e r deta il in the articles by 
iBoussol (120071) : iPadmanabhaJ (|2007l) in this special vol- 
ume. 

To paraphrase Winston Churchill speaking about 
democracy, it may be that the cosmological constant is 
the worst form of accelerating physics, except for all 
those other forms that have been tried from time to 
time. Nevertheless, this article addresses those other 
forms, specifically dynamical physics that aims to ame- 
liorate the coincidence, and/or fine tuning, problems. 
We concentrate on the dynamics, the time evolution 
of the cosmological expansion physics, (mostly) from 
a c anonical scalar field, given th e nam e quintessence. 
See ICopeland. Sami. fc Tsuiikawal (2006) for a particle 
physics perspective. 

Section HT1 provides a brief historical perspective on the 
development of quintessence theories. Section lnTl rcvicws 
key elements of the dynamics of quintessence and the 
physical origins of structure in the phase space, defin- 
ing classes of models. Efficient representation of the dy- 
namical behavior through parameterization or principal 



component analysis is discussed in Section IPVl and we 
investigate in detail thawing models, those which ap- 
proach cosmological constant behavior, in Section IVIII 
In Section [Vj we consider a selection of dynamical mod- 
els beyond standard quintessence, and briefly mention 
the effects of expansion dynamics on growth of cosmic 
structure in Section IVT1 We conclude in Section [Villi 



II. ORIGINS OF QUINTESSENCE 

The role of a dynamical scalar field for recent acceler- 
ation of the cosmic expansion certainly owes a debt to 
the use of rolling scalar fields for early universe inflation. 
A scalar field, and more generally a negative equation 
of state, were implemented as a substitute for the cos- 
mological constant in a flurry of ac tivity in the 1980s. 
On the theoretical side. iLindg ( 19871 ) proposed a simple 
extension from the flat potential of the cosmological con- 
stant to a tilted, linear potential, that releases the field 
to roll when the expansion rate of the universe decreases 
sufficiently, what is now called a t hawing field. In 1988, 
two nearly simultaneou s papers by I Wetterichl ( 19881 ) and 
iRatra fc Peebl es ( 1988) described in more detail cosmol- 
ogy in the presence of a quintessence field. 

At the same time, considerable work on the phe- 
nomenology of energy density components with an ar- 
bitrary (including negative) pressure to density, or equa - 
tion of state, ratio was being carried out. IWagoner (11986 ) 
discussed such generalized cosmolog y, and lLinderl ljl988af) 
then followed up on this with detailed investigation of a 
variety of cosmological probes of additional components 
with arbitrary equation of state. These included tests 
of the expansion dynamics through distance, age, vol- 
ume, and abundance measurements. Particular attention 
was paid to light propagation in such a generalized cos- 
mology, including poss ible inhomogeneities in the com- 
ponents (|Linder[[l988oT ) 



some re s ults o ccurred earlier in 
the unpublished thesis of iKavserl (|l985l )). General equa- 
tions of state had been considered in a formal way for the 
growth of structur e with in linear perturbation theory by 
Kodama fc Sasakil (|1984l ) . Implications of ge nera l equa - 
tion s of state for growth were presented in iFrvi ( 1985f ) 
and|Lmder| (|l988dY 
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Thus high energy physics theory and cosmology were 
all ready in the 1980s for data exploring the expan- 
sion and growth histories of th e universe. It took an- 
other 10 years fo r observations (jPerlmutter et all Il999t 
iRiess et all Il998f) to make the astonishing breakthrough 
that turned these speculations into a central subject of 
research into our understanding of gravitation, quantum 
physics, cosmology, and the fate of the universe. 



III. THE QUINTESSENCE OF DYNAMICS 
A. Scalar Field Basics 

If we view the cosmological constant as a quantum 
zeropoint energy corresponding to the ground state of 
harmonic modes of a field filling space, we can picture 
this as an array of identical springs, motionless and each 
stretched to the same length. By contrast, a scalar field 
would be a dynamical version of this, with the springs os- 
cillating in time and having different lengths at different 
points in space. That is, a scalar field is a very simple 
quantity, a magnitude at each point in space. One can 
literally picture it as a field: a field of grass where each 
stalk may have been mown to a different height (a vector 
field could then be a field of trampled grass, where each 
stalk has a length and a direction in which it lies). 

For quintessence, we take a scalar field cj) minimally 
coupled, i.e. feeling only gravity, passively through the 
spacetime curvature, and a self-interaction described by 
the scalar field potential V(4>). Moreover, we consider the 
kinetic contribution to the Lagrangian (the "bouncing of 
the springs" ) to be canonical, i.e. involving only a term 
linear in the kinetic energy of the field. (We briefly dis- 
cuss relaxing these conditions in SjV]) So the Lagrangian 
is about as simple as possible: 



C^ = ^d^ct>-V{<t>). 



(1) 



Through the Noether prescription we define an energy- 
momentum tensor 



2 6(y/=gC) 
/=g Sg^ 



(2) 



where t^,„ is the metric and g its determinant. Compar- 
ing the result for a homogeneous and isotropic spacetime 
to the perfect fluid form allows identification of the en- 
ergy density and pressure: 



(3) 
(4) 



Because late time acceleration requires a very light scalar 
field, with effective mass of order the Hubble parameter, 
the Compton wavelength of the field will be of order or 
larger than the Hubble scale and so the field is expected 



to be spatially smooth within the Hubble scale. There- 
fore we neglect the spatial gradient terms in the energy 
density and pressure. These quantities can be put into 
the usual Friedmann equations to solve for the expansion 
history of the scale factor vs. time, ait), from the Hubble 
parameter H = a/a and acceleration d. 

Because both the energy density and pressure enter the 
equations, it is convenient to define an equation of state 
ratio, 



(5) 



which is generally time varying. When we refer to dy- 
namical fields, we generally mean time- varying w, i.e. 
W ^ constant. (Although the energy density of constant 
w models varies with time, this happens as well with 
matter or a frozen network of cosmic strings, say, and so 
does not capture the flavor of "dynamics".) 

The equation of motion for the scalar field is the Klein- 
Gordon equation 



3H</> = -dV/dcj), 



(6) 



and is interchangeable with the continuity equation. For 
example, multiplying through by <f> gives the sequence 



[0 2 /2]-+6 J ff[^ 2 /2] = -V 



p<j> - V + 3H(p<t, + p^) 

4f^- = -3(P0 +P<a) 
din a 



-V 

-3,00 (1 + w) 



(7) 
(8) 

(9) 



where we have turned Eqs. ©-(HI) around to use 

V = [p (j> -p <j> )/2 = P<i> {l-w)/2 (10) 
K = 4> 2 /2 = ( P4 ,+ P<t> )/2 = p^{l + w)/2. (11) 

From the above equations we can formally go back and 
forth from the field description to the fluid description or 
equation of state. From Eqs. ©-{J) we see that 



K-V 
K + V' 



(12) 



so for some specified theory we can calculate the equation 
of state and then the effects on the cosmological expan- 
sion. The other direction, starting from observations of 
the cosmological expansion, is slightly more complicated: 



p^(a) — Cl w p c exp -^3 J dhxa [1 + w(a)] j (13) 

dlnaH- 1 ^ P<j) {a) [1 + w(a)} (14) 
P<t> {a) [1 - w(a)}/2 (15) 



0(a) 

V(a) 
K{a) 



<f> 2 /2 = P<t> {a) [l + w(a)]/2. 



(16) 



Such reconstruction of the scalar field physics is made dif- 
ficult by a number of issues: noisiness of measurements 
of the expansion, translation from the measured quan- 
tity to density or equation of state through one or two 
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derivatives, and finite range of scale factor, or redshift 
z = a^ 1 — 1, coverage. In particular, from the last of the 
equations above we see that 

4> = Ml + ™)] 1/2 £ HM P (1 + w) 1 / 2 , (17) 

so for cases when l+w <C 1 (as seems to be implied by ob- 
servations), only a small region of the scalar field physics, 
Acf> ~ 4>/H <C Mp, can be probed. All these issues to- 
gether makes reconstruction problematic, and we do not 
co nsider it further. (For attem pts to carry it through, 
see lSahni fc Starobinskyl ( 20061 ) and references therein.) 

While we cannot reconstruct in detail the scalar field 
potential, we can derive considerable insight into the ac- 
celerating physics from study of its dynamics. We can 
guess from the spring picture at the beginning of this 
section that there will be at least three basic quanti- 
ties we want to know: how much energy is there in 
the field, how springy is it, and how stretchy are the 
springs? The energy density is conveniently written 
in terms of the dimensionless density Q w — p^j p c , where 
p c = 3Hq I (8irG) is the critical density. For a spatially 
flat universe, fl w = 1 — il m , where f2 m is the dimen- 
sionless matter density. The analog of the springiness 
is how spacetime curvature reacts to the accelerating 
component; the passive gravitational mass is given by 
p + 3p, with acceleration induced by a component pos- 
sessing p < — p/3, or w < —1/3. So we can regard w as a 
measure of the springiness. As the universe expands, the 
springs react, changing their springiness, like stretching 
the coils of a spring. This time variation can be taken as 
w' = dw/dhxa = w/H. Thus we are primarily interested 
in Q w , w, w' . The last two quantities give a phase space 
for the dynamics which we will see is enlightening. 



B. General Dynamical Behavior 

Scalar fields can at any epoch have one of four behav- 
iors. Their rolling can be fast, slow, more or less steady, 
or oscillatory. 

Fast roll: Fast rollers have kinetic energy exceeding 
their potential energy, and so by Eq. (fT2")) have w > 0. 
These clearly do not act to accelerate the cosmic expan- 
sion, but a fast roll epoch ( "kination" ) is a characteristic 
of tracker models, which follow attractor trajectories in 
their dynamics such that at certain epochs their equation 
of state is determined by the dominant energy density 
component of the universe. Because of the fast roll, the 
scalar field can rapidly decrease its energy density from 
an initial, early universe value near the "natural" Planck 
scale to a much smaller value that will make it suitable 
for the observed present energy density. Due to the at- 
tractor solution for the dynamics, for certain forms of 
the potential, there is a large variety of initial conditions 
- "basin of attraction" - that can deliver a reasonable 
present energy density, thus addressing the fine tuning 
problem of the cosmological constant. Of course the field 
must leave both the fast roll regime and the tracking 



regime if it is to cause acceleration and dominate the en- 
ergy density, so the coincidence problem is not completely 
solved. In particular, tracking fields have difficulty reach- 
ing equations of state w ^5 —0.7, in tension with obser- 
vations, and so are no longer considered front runners 
for explaining the acceleration . For more on tracker s 
(and the earlier "t r acers" ), see iFerreira fc Jovcd (11997T): 
iLiddle fc Scherrerl Jl999h: ISteinhardt. Wang, fc Zlatevl 



(jl999l ); IZlatev, Wang, fc Steinhardtl (|l999t ) 



Slow roll: When the kinetic energy is much smaller 
than the potential energy, the equation of state is 
strongly negative, w w -1. Of course this only leads to 
acceleration of the expansion if the dark energy also dom- 
inates the energy density. The field is nearly frozen, and 
the dark energy density is nearly constant (while matter 
and radiation are rapidly diluting due to the expansion) , 
so it would eventually come to dominate the universe 
if nothing else changed. Note that because matter is 
not negligible, even today, a field we think of as slowly 
rolling, w ~ — 1, may well n ot have a small value for 
V'/V (see, e.g. lLinderl (12006^ )1. which is a conventional 
slow roll parameter for inflation (where the accelerating 
component is completely dominant). Quintessence mod- 
els that always have the potential dominating over the 
kinetic term encounter the same fine tuning and coinci- 
dence problems as the cosmological constant, lacking the 
basin of attraction of tracker models. Thus generically, 
we want a combination of fast and slow roll behavior for 
a successful model. 

Steady roll: Refer ring to the original quintessence 
model of lLindd (|l987h using a linear potential, this cat- 
egory is somewhat of a misnomer since the field does 
have fast and slow roll epochs over its entire history. 
However, the linear potential model does have a con- 
stant right hand side of the Klein-Gordon equation of 
motion, and for a long time the dynamics stays reason- 
ably close to the line where the field acceleration <p (not 
the cosmic acceleration a) is zero (see i jlll.Cl below). This 
model is not only the simplest generalization of the cos- 
mological constant but is also interesting in its overall 
history. It starts generically from a frozen, cosmological 
constant-like state due to Hubble friction, then thaws and 
rolls down the potential. However, because the potential 
has no minimum, the field rolls into territory where the 
potential goes negative, which actually leads to a col- 
lapsing universe, rather than an accelerating expansion. 
These models ther efore have a fin i te future history, with a 
"doo msday time" ( Kallosh et al.l . [20031 iKratochvil et all 
12004ft . 

Oscillation: Common potentials in renormalizablc 
field theories include V((f>) ~ <j) n , which have a minimum 
for n even. While the field will have a conventional rolling 
stage, eventually it will reach the minimum and oscillate 
around it. If the period for oscillation is much smaller 
than the Hubble time (as is general ly the case) th en the 
effective equation of state becomes ( Turnerl . 11983) 



w 



(18) 
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For a quadratic potential, the field acts like nonrelativis- 
tic matter, and for a quartic potential it acts like radia- 
tion. 

One intriguing example of such a field is the ax- 
ion, or more generally pseudo-Nambu Goldstone bosons 
(PNGB). If we consider them during the regime when 
they are still rolling rather than oscillating, they can 
accelerate the expansion, though this acceleration will 
eventually fade aw ay as the field evolves to its oscillatory, 
matter-like phase ( Frieman et all [l995l ). PNGB poten- 
tials are also radiatively stable against quantum correc- 
tions, unlike an ad hoc V(cj>) that might be written down 
but then acquire a non-zero ground state (cosmological 
constant) and distortion of its shape. Thus the physics 
of such pseudoscalar fields offers some promise for a fun- 
damental, high energy physics origin rather than merely 
a low energy effective potential. The PNGB potential 
looks like 



V(4>) = V [l + cos(<}>/f)], 



(19) 



where / is a symmetry energy scale. Because the 
potential is nonmonotonic and the slope of the po- 
tential changes from concave to convex, a num- 
ber of interesting effects can arise, such as mimick- 
ing super- negative e quations of state w < — 1 and 
nontrivial dynamics (fCsaki. Kaloper. fc Terning] . 20061 : 



lHall. Nomura, fc OliverL 120051 : iKaloper fc SorboL l2006h . 
For a complex field, one has degrees of freedom in both 
the modulus and the phase, and researchers have con- 
sidered making on e act as dark energy and the other as 
dark matter (e.g. iMainini fc Bonomettol . |2004( ). or one 
giving recent acceleration a n d one early universe inflation 
(e.g. iRosenfeld fc Frieman . 2005h. Other elaborations 



include spintessence ( Boyle. Caldwell, fc Kamionkowskil 
l2002t iGu fc Hwanel . l200lh . 



C. Fundamental Modes of Dynamics 

While in the previous subsection we considered the be- 
havior of the scalar field dynamics at any one moment, 
considerably more insight comes from investigating the 
overall dynamical history given by the trajectory through 
phase space. In particular, we will be interested not only 
in the present characteristics, but the asymptotic past 
and future states. 

By examining the physical impact of the three different 
terms in the Klein-Gordon equation ([6]) we can identify 
boundaries in the phase space corresponding to different 
physical conditions. 

• Phantom line: This line separ ates physics obey- 
ing t he null energy condition (jHawking fc Ellis] . 
1973), p + p > (w > —1), from physics violat- 
ing it. Also, consider the friction term 3H<fi. From 
Eq. (jlip one sees that where the sign of this term 
changes, i.e. <j> — as the field stops rolling in one 



direction (and possibly begins rolling in another), 
corresponds to 



1. 



(20) 



Canonically the field has w > — 1 but there are 
various mechanisms (see for achieving w < 
— 1, what is re ferred to as the phantom regime 
( Caldwell [20021) . 



• Null line: Consider the forcing term of the poten- 
tial slope. When the field rolls down the potential, 
V < 0, this corresponds to 



>/ > -3(1 -w 2 ), 



(21) 



where we have used Eqs. ©-(HU to convert the 
variables V and <f> to w, w' . If the field has a (non- 
canical) negative kinetic energy so it rolls up the 
potential then the inequality flips but at the same 
time the sign of w changes so ui < —1 (one can 
think of this as the energy density increasing with 
time, following Eq. [5]). Thus the null line passes 
smoothly through the point [w, w') = (—1,0). 

• Coasting line: Consider the acceleration term <j>. 
Generally, at late times, the field accelerates due 
to the potential forcing dominating over the fric- 
tion, or decelerates if the friction dominates over 
the potential slope (note this should not be con- 
fused with the acceleration of the cosmic expan- 
sion, which holds in either case if w is sufficiently 
negative). Again from Eq. (TTTjl the dividing line 
between these dynamics, where the field is freely 
coasting at constant velocity 4>, is 



w' = 3(1 + wf 



(22) 



with w' greater (smaller) than this for field accel- 
eration (deceleration) . 

These three boundaries give general physical divisions 
for the dynamical behavior of the field. The general equa- 
tion relating the phase space variables can be derived by 
taking the derivative of Eq. (|10p and using the continuity 
equation © to obtain 



-3(1 - w 2 ) - (1 - w){\ + w) 1/2 \j 



8tt 



V 



(23) 

We can readily verify the null line corresponds to = 
(and one can specialize to the coasting line with 
a little more effort). The se conditions were defined 
in ICaldwell fc Linderl (120051 ) and d eveloped further in 
IScherreri (|2006h and lLindeiT(l2006bl ). The last reference 
in particular goes into more detail about the derivation 
and the effect of the ratios of different terms in the Klein- 
Gordon equation, as well as "slow roll" parameters of the 
potential. 

Without the need for quantitative analysis of the ra- 
tios of Klein-Gordon terms, one can broadly understand 
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the dynamics by examining the relati ve dominance of the 
drivin g vs. dragging terms, following ICaldwell fc Linderl 
[f the Hubble friction dominates at early times, 
then the field will be pinned and act like a static cos- 
mological constant. As the cosmic expansion reduces the 
Hubble parameter, eventually the potential slope induces 
the field to begin rolling: such models are said to be 
thawing, and their dynamics in phase space shows them 
"leaving A", moving to less negative w with positive w' . 
In particular, fields that thaw during the matter domi- 
nated epoch leave A along the track w' = 3(1 4- w). 1 As 
the matter domination wanes, the trajectory will curve 
according to the driving force from the potential slope; 
since the potential (eventually) becomes less steep as it 
approaches the minimum, the field acceleration decreases 
and the curve is toward the coasting line, i.e. smaller w' . 
For broad classes of potentials the condition that dark en- 
ergy not completely dominate the energy density of the 
universe by the present means that thawing fields are still 
accelerating along the potential and the dynamics has a 
lower bound roughly given by w' > 1 + w (for fl w < 0.8 
and w < —0.8). Thus the thawing region of phase space 
is defined by a dynamical history 



history 



1 



< 



/ < 3(1 +w). 



(24) 



The alternative is that the potential forcing dominates 
over the Hubble drag at early times, i.e. the potential 
is sufficiently steep to overcome the friction from cosmic 
expansion. Such fields will look different from the cos- 
mological constant at early times. Certain forms of po- 
tential possess special attractor properties, as discussed 
in the previous subsection, that during the matter dom- 
inated epoch cause the scalar field dynamics to have a 
constant equation of state determined by the background 
expansion. As the dark energy density becomes relatively 
more important, these fields will depart from their track- 
ing behavior and roll according to the dynamics of their 
potential. As the field rolls toward the minimum, de- 
celerating in its motion (lying below the coasting line), 
gradually approaching asymptotically a static cosmolog- 
ical constant state, it is said to be freezing. In its "ap- 
proaching A" , the field c o ntribut es an energy density 
p w ~ H 2< - 1+w \ but iLinderl (|2006bl) showed that any H a 
model approaching w = — 1 does so along the asymptotic 
trajectory w' = 3w(l + w). Conversely, since dark en- 
ergy dominates (though not fully) today, the field must 
have departed its matter dominated tracking behavior 
and moved some distance away from the constant w line. 
For broad classes of potential this leads to a present value 
w' < 0.2w(l + w) (for Q w > 0.6 and w < -0.8). Thus the 
freezing region of phase space is defined by a dynamical 



1 Fields whose initial conditions <pi are fine tuned can avoid this. 
Also, if the potential driving term is very large, for example in 
PNGB fields with symmetry energy scale / <C Mp, then one can 
have vi' > 3(1 + w). 



< w' < 3w(l + w), 



(25) 



with the present value of w' more tightly restricted. 

Figure Q] illustrates the three critical dividing lines of 
the phantom, null, and coasting curves in the dynam- 
ical phase space. In addition it shows the upper and 
lower boundaries of the thawing and freezing regions. 
Note that the lower boundary of the freezing region coin- 
cides with the constant pressure curve (with an adiabatic 
sound speed c\ = 0) discussed in SjV] 
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FIG. 1 The dynamical phase space w-w' is divided by three 
curves defined by physical conditions: the phantom line 
w — —1, the null line w' = —3(1 — w 2 ) following from a flat 
potential, and the coasting line w' — 3(1 + w) 2 following from 
constant field velocity. These extend across the phase space. 
In addition, canonical dynamics leads to the distinct regions 
of the thawing regime bounded by the red dotted lines and the 
freezing regime bounded between the green dot-dashed curve 
and the blue dashed curve (the latter given by the constant 
pressure condition). 



Comparing Eqs. (|24|) and (|25|) . we see that they define 
narrow, distinct regions in the phase space where scalar 
field theories obeying a combination of theoretical and 
observational conditions lie. In particular, there are fairly 
strongly physically motivated outer boundaries defining 
the extremes of w' . The exact inner boundaries are more 
a function of empirical constraints on the present ex- 
pansion, but there is a distinct intermediary zone un- 
favorable for habitation. This "desert" lies around the 
coasting line: only highly fine tuned models would, af- 
ter the many e-folds of cosmic expansion influencing the 
scalar field equation of motion, find themselves almost 
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perfectly balanced between field acceleration and decel- 
eration, (p ps 0. 

Two important implications of the physical division 
into distinct thawing and freezing regions are for the 
questions of observationally distinguishing dynamical 
dark energy from A and distinguishing the physical ori- 
gin of the dark energy (e.g. field theories with thawing 
vs. freezing characteristics). Because of the degeneracy 
directions of essentially all co smol ogical p r obes (see the 
articles bv lLeibundeull (|2007ri and lNicholsl (|2007h in this 
volume), the entire thawing region is difficult to distin- 
guish from the cosmological constant if the data is only 
at the sensitivity level of a constant, or time averaged, 
w. For example, the entire thawing region would give 
an apparent (w) ~ —1 ± 0.05. Thus experiments sen- 
sitive to w' are necessary for deciding between this half 
of the dynamical phase space and the cosmological con- 
stant. For distinguishing between the classes of effective 
field theories, one would like to have cosmological sensi- 
tivity to the time variation of a(w') < 2(1 + w) to resolve 
the separation between the thawing and freezing regions. 
For indepth discussion of mappi ng the cosmic expansion 
history, see the review article bv lLinderl ((20071 ) . 



D. More Complicated Dynamics 

In the previous subsection we gave physical motiva- 
tions for bounded regions in phase space but we empha- 
size these are based on a combination of generic behavior 
and empirical data, not an absolute exclusion of other 
possible behaviors. In particular, they relied on a stan- 
dard matter dominated epoch at high redshift, canonical 
scalar fields, avoidance of fine tuned initial conditions and 
potential shapes, and "fundamental modes" of dynamics. 
We discuss extension of the dynamics to beyond canoni- 
cal scalar fields in SjVj here we consider initial conditions 
and fundamental modes. 

Initial conditions on the scalar field dynamics are quite 
important, e.g. one could consider a field so perfectly bal- 
anced on a maximum of its potential that it only starts 
rolling yesterday, or a field that has recently passed a 
minimum of its potential and is now rolling uphill, or a 
field with kinetic and potential energies exactly crafted 
so the dynamics is missing (constant equation of state) 
or is coasting. Physics does not forbid any of these a 
priori, but our sense of naturalness disfavors them. If 
dynamical conditions are set by hand at recent times, 
rather than the field settling into an evolution follow- 
ing its equation of motion over many e-folds in the early 
universe and then a matter dominate d epoch, then vir- 
tually arbitrary behavior can r esult ( Huterer &: Peirisl 
l2007t iLi. Holz. fc Cooravi l2007h . One could fine tune 
the field such that one does not extract general physi- 
cal precepts on the dynamics, but rather the phase space 
trajectories would spell out your name. 

Under the physics of field evolution through the cosmic 
expansion history, including a matter dominated epoch, 



the phase space structure described in the previous sub- 
section generically holds. One further necessary ingredi- 
ent is that we are talking about fundamental modes, or 
"atoms" , of the dark energy - the quintessence of dynam- 
ics. If one combines multiple elements together, such as a 
scalar field plus a cosmological constant, or plus matter, 
or plus another scalar field, then one can indeed break 
the physical boundaries (just as multifield inflation can 
break consistency relations and other basic predictions) . 
That is, the phase space structure applies to the dynam- 
ics of a single, fundamental field, not an effective field of 
multiple origins. 

We can investigate this further by examining the ef- 
fect on the equation of state when multiple elements 
are combined. For the simplest approach, we consider 
adding together two components: a canonical scalar field 
plus either a cosmological constant, a matter component 
(e.g. misestimation of Q m or dark energy contribution to 
dark matter), or another scalar field. 

The effect of combining two such noninteracting com- 
ponents is given by an effective dynamical equation of 
state 



W c g = Wi 



SH 2 



5H( 



SH 2 



U>2 



SH 2 



5H 2 



SHI 



(26) 



where SHf is the contribution of component i to the 
Fricdmann equation. This approach was used to first 
point out phanto m crossing, ev olution a cross w = —1, by 
two scalar fields (|Linderl . [2003 ) (also see[Hu| (|2005l )). The 
dynamics is affected as 



3w e ff(l 



Wefi) 



SH 2 



3w;i(l + wi)] 



:tt2 

— 2V 2 -3«*(1 +«*)], (27) 



where S = SH 2 + 5H 2 . Note that two constant pressure 
components (where w[ = 3z«i(l + Wi)) add without af- 
fecting the dynamics. In particular, any combination of 
matter plus A keeps the same trajectory, just moving the 
position along the track. 

Furthermore, this formula implies that the sum of com- 
ponents, each of which lies on the same side of the curve 
w' = 3w(l + w), has effective dynamics doing likewise. 
For example, two kinetic k-essence components give an 
effective dynamics that is still kinetic k-essence-like. Sim- 
ilarly, the null condition w' > — 3(1 — w 2 ) cannot be over- 
come by summing components obeying w[ > —3(1 — wf). 
Other than respecting these two boundaries, the dynam- 
ics can change significantly on combining components. 

To an initial thawing scalar field we add either a cos- 
mological constant component, a matter component, or 
a freezing field. Figure shows that such combinations, 
as opposed to the fundamental modes or "atoms" we 
discussed in the previous subsection, do not adhere to 
the restricted thawing and freezing regions of the phase 
space. Convolutions of different physics can drastically 
differ from those fundamental behaviors. 
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FIG. 2 Dynamics involving combination of physics can violate 
the fundamental phase space regions. To the original thawing 
scalar field trajectory (solid black), we add a cosmological 
constant (+A), extraneous matter or quartessence component 
(+m), or freezing scalar field (+V). We fix wq = —0.8 for the 
fields and take the total dimensionless dark energy density 
to be 0.7. For the second component of A or V we take 
^2 = 0.1 (darker, black) or 0.35 (lighter, red); for included 
matter fl+ m = 0.01. Curve endpoints correspond to z = 0, 
with x's at z — 1. 



Adding a freezing field to a thawing field dramatically 
alters the trajectory, since at early times the freezing field 
will dominate. (Adding extra components to a dominant 
freezer has less effect.) The phase space tracks therefore 
start off in the freezing regime but curve up toward the 
thawing regime, possibly lying today in the desert re- 
gion between the two regimes. A cosmological constant 
rotates the dynamics to ward w' = and d r aws it in to- 
ward w = -1 (see also ICaldwell fe Linderl 1)20051 )): this 
does not generally move a thawing field out of the thaw- 
ing region. Including a matter like component with the 
thawing field has the most severe effect. Adding a mere 
0.01 in dimensionless matter-like energy density alters 
the track wildly - this points up strongly the dangers 
in attempted direct reconstruction of the dynamics from 
H(z) or the distance-redshift relation. Misestimation of 
Q m by 0.01 will completely distort the true dark energy 
dynamics. 



IV. DESCRIBING THE DYNAMICS 

The phase space dynamics discussed in the previous 
section presents the dark energy physics in terms of 



a function w(a) and its derivative it/, describing the 
"springiness" and "stretchiness" of the spacctimc in reac- 
tion to the dark energy. Each theoretical model presents 
its particular description of the function and we can check 
each against the data to determine whether the model 
fits. However, there are 10 x theoretical forms (potentials 
or equation of state functions) already postulated, each 
with their own parameters. Moreover, we would like to 
predict the results of experiments, or design experiments, 
more generally than for a given theory or set of existing 
theories. 

This shows the need for a model independent ap- 
proach, based on a parametrization of the equation of 
state function or a similar quantity. Because we want the 
parametrization to stay close to the underlying physics, 
of which both the dark energy density and pressure 
enter, we concentrate on the pressure to density ra- 
tio, or equation of state ratio. However parametriza- 
tion of other quantities such as distances, Hubble pa- 
ram eter, or density alone hav e been considered (see, 
e.g.. lSahni fc Starobinskvl ( 20061 ) and references therein); 
two cautions should be stated about this route: cer- 
tain for ms bias the ex t ractio n of the underly ing physics, 
see e.g. iJonsson et al. ] (I2004D : iLinderl (|2006d) . and if one 
eventually wants the equation of state then one is forced 
to take numerical derivatives of a quantity extracted from 
noisy data. 

Numerous parametrizations exist for the equation of 
state w(a) but the vast majority are purely ad hoc. 
We here consider a very few that are phenomenologi- 
cal in the best sense, i.e. generalized from the behavior 
of physically motivated sets of models. From the previ- 
ous section we have seen that a single parameter model, 
i.e. w — constant, involves highly fine tuned physics to 
remove the dynamics. While one way out of this is to 
invoke a physical symmetry, such as a topological defect 
origin, which can produce w — —N/3 for a frozen net- 
work of A^-dimensional defects ( e .g. N = 2 domain walls 
(IZel'dovich. Kobzarev. fc Okunl . I1975D or N = 1 light 
cosmic strings ( Vilenkinl . 1 1 9 84 ) ) , such values are not con- 
sistent with data. 

This leads us to two parameter models as the next 
simplest alternative. The parametrization 



w(a) = wo + w a (l - a), 



(28) 



where wq is the value today (a = 1) and w a is a measure 
of the time variation w' , is widely used in the literature. 
It is important to realize that it is in no way a mathe- 
matical expansion a bout the present: neither its impor- 
tant introduction bv lChevallier fc Polarski (l200ll) nor the 
physical foundation work by iLinder (2003a) employed 
a Taylor expansion, nor would that be mathematically 
convergent. Therefore w a is not an expansion parame- 
ter about z = 0, but rather a fit parameter describing 
the ove r all tim e variation w' . The original convention 
([Linderl . l2003al ) giving the best description is 



w a = (-u//a)| 2 



-2w'(z = 1). 



(29) 
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iLind er ( 2003aHJ) give several physical supports for the 
w a parametrization: 1) excellent approximation to the 
exact field equations for a broad range of fundamental 
or straightforward scalar field potentials, 2) well behaved 
at both low and high redshift, 3) robust against bias, 
e.g. if one extends the form to further parameters, the 
wo, w a parameter values estimated are not strongly af- 
fected, 4) model independence. For example, a SUGRA 
inspired model that evolves from w(a < 1) « —0.2 to 
w = —0.82 - a substantial variation - has its equation 
of state reproduced to within 3% back to z — 1.7 and 
the distance-redshift relation in such a cosmology is ac- 
curately matched to 0.2% back to CMB last scattering 
by w = -0.82, w a = 0.58. 

Of course a two parameter description cannot describe 
all possible dynamics; in particular it begins to break 
down for rapid transitions in the equation of state or 
oscillations. However, for the fundamental modes high- 
lighted in the previous section it serves as an excellent, 
broad (i.e. model independent, good for both thawing 
and freezing) parametrization of the physically favored 
dynamics. 

Another two parameter form, which is motivated from 
the energy density rather than the equation of s tate, is 
the bending parametrization of IWetterichl ( 2004h . This 
was designed to describe early dark energy models where 
at high redshift (near the CMB last scattering surface, 
z w 10 3 ) the scalar field component has nonnegligible 
energy density (though it is then acting in a decelerating, 
rather than accelerating, manner on the expansion, so it 
is not exactly dark energy). The bending form has 



In 



Q m (a) 
w(a) 



— Ro — 



3wq In a 
1 — b In a 
w 



(l-&lna) 2 ' 



(30) 
(31) 



where Ro = ln^" 1 — 1) and b is related to the early dark 
energy density. The dynamics of this parametrization is 
that in the past it approaches w = 0, w' = (i.e. a finite 
dark energy density that acts like matter) , at some future 



time a* = e 1 ^ it runs to w — 



= — oo, and then 



returns along the same trajectory to w — 0, w' = in 
the further future. The phase space track is defined by 
w' = 2bw (w/w ) 3 / 2 - At any given time in the past the 
variation must be slower than w' = — (8/27)z«o/ lna. 

A generalization of the w a form to three parame - 
ters was put forward bv lRapetti. Allen, fe Wellerl (|2005h . 
This eases the property of the w a form where the param- 
eter w a plays two roles: it describes the characteristic 
time variation w 1 but it also determines the asymptotic 
past value of w(a <C 1) — > wo + w a - The extended form 
has 



w(a) 



w p z + w z t 
z + z t 



(32) 



where w p is the asymptotic past value and z t is the tran- 
sition redshift. When z t — I, this reduces to the w a 



parametrization. The phase space dynamics is a parabola 
from (w, w') — (w p ,0) to (u>/,0), crossing w = — 1 if 
w < w p . 

To describe a monotonic w(a) which transitions 
smoothly from some asymptotic past value w p to some 
asymptotic future value Wt requires a minimum of four 
parameters: w p , Wf, the epoch of transition et t , and a 
rapidity parameter r. (Note that the previous models 
are not bounded in the future; this is not overly worri- 
some because we have no data on the expansion future.) 
Such forms are particularly successful in describing track- 
ing models which have both asymptotic past and future 
equations of state. The transition can be described by 
many functional forms, but the two most common four 
parameter equations of state b oth adopt "Fermi-Dirac" 
trans itions. The kink model ( Corasaniti fe Copelandl . 
2003) takes this in scale factor a, obtaining 



w(a) = w + (w m - w Q ) 



1 _|_ e o«/A i _ e (l-o)/A 
1 _ e l/A i _|_ e (o t -o)/A : 



(33) 



where w m is the asymptotic value in the matter domi- 
nated era a nd A is related to the r apidity, while the e- 
fold model (jLinder fc Hutererl . 120051 ) does the transition 
in the expansion e-fold factor In a, obtaining 



w(a) — Wf + 



W p — Wf 

1 + (a/a t y/^ 



(34) 



One of the advantages of the e-fold model is that it allows 
an analytic expression for the Hubble parameter H(a). 

One could continue developing more complicated forms 
but sadly even the next generation of experiments will 
not be able to constrain s tringently more than two equa- 
tion of state parameters ( Linder fc Hutererl . [2005) . This 
conclusion holds whether dealing with parameters per se 
or principal components (see below). Happily, the w a 
parametrization is quite satisfactory in giving a model 
independent, good approximation to the dynamics. 

Nevertheless, let us briefly consider principal com- 
ponent analysis (PCA). This approach attempts to 
gain some independence from the particular form of 
parametrization, letting the data define the best con- 
strained combi nation of information . Th is is a valu- 
able tool; see iHuterer fc Starkmanl ( 2003h for its de- 
velopment for the dark energy equation of state, and 
IHuterer fc Cooravl (l2005h for an adaptation localizing the 
principal components in redshift. PCA has the advantage 
over parametric forms in its nonparametric flavor, and in 
specifying what a particular survey measures best, how- 
ever its results are dependent on ingredients other than 
the underlying physics: the type of cosmological probe, 
the details of the data, the fiducial cosmology, and priors. 
That is, a principal component derived from one specific 
experiment is not exactly comparable to a principal com- 
ponent from another experiment, or the same experiment 
over a different redshift range. By contrast, wq and w a , 
say, mean the same thing regardless of probe, survey, cos- 
mology, or priors. (We are talking about the meaning of 



9 



the variables, not the estimation of the fit values.) Thus, 
PCA is likely to be of most use as a complementary tool 
alongside parametric fits. 

Note there has been some confusion in the literature 
regarding the accuracy of PCA fits, with some claims that 
more than two principal components can be stringently 
fit by next generation experiments. In the analyses where 
there appear to be more than two well fit parameters, 
this arises from consideration only of low noise in the 
component coefficients ccj, e.g. a{o.i), not high signal to 
noise criteria cr(aj)/aj. 

So we appear restricted to two parameters for our equa- 
tion of state description. However, a tilt from the cos- 
mological constant value, 1 + w, and a time variation, 
w', contain rich information on the physics responsible 
for the acceleration of the universe. Given we have only 
two parameters, are we sure that wo and w a represent 
the best, model independent parameters? No, we have 
no guarantee of this and we should continually be on the 
lookout for improvements, though to date wq, w a have 
served extremely well. 

One idea for an alternate parameter involves the so- 
called pivot or minimum variance equation of state w p . 
This is the equation of state at the scale factor a p where 
the variance o 2 (w(a)) is minimized, i.e. w p = w(a p ). 
Note that w p is also decorrelated with w a , with zero co- 
variance between their estimations, but this holds only 
due to the specific linear dependence of the equation of 
state w(a) on w a ] generally the minimum variance value 
is not decorrelated with other equation of state param- 
eters. The pivot parameter possesses many of the same 
issues as the PCA approach: lack of an invariant physi- 
cal meaning due to dependence on probe, survey, model, 
and priors. It is sometimes useful however for the narrow 
question of whether the data are consistent with a cos- 
mological constant cosmology (in one direction, at least; 
one can find w p = — 1 yet have dynamical dark energy). 
Note for thawing models the deviation 1 + w is greatest 
at z = so a par ameter a t z v m ay not be optimal even 
for this question. iLinderl ( 2006d) showed that generally 



w p is more subject to bias than either wo or w a . 

Another suggestion for alternate parametrization 
i nvolves either s o-called statefinder variables (r, s) 
(|Sahni et al.l . l2003h or combinations of derivatives of the 
cosmic scale factor such as the decele ration parameter 
q = —ad/ a 2 and jerk j = a 2 'a/a 3 ( Blandford et all 
120051 ). Note that either parametrization convolves the 
equation of state parameters with the energy density: 



1 3 / \ 



(35) 



1 - -il w (a) K - 3w(l + to)] =q + 2q 2 - q'{36) 



(r is the same as j, and s = [3w(l + w) — w']/(3w) = 
c 2 (l + 11)) /w, where c 2 is the adiabatic sound speed.) 
These approaches also conflate different physics: j — 1, 
for example, corresponds to an Einstein-de Sitter pure 
matter universe, or a de Sitter pure cosmological con- 



stant universe, or any model that instantaneously lies on 
the w' — 3w(l + w) line. Of course interpreting q and 
j as a Taylor expansion about the present expansion be- 
havior would restrict their usage to z < 1. Also note 
that while the scale factor can be viewed as a kinemat- 
ical quantity (e.g. no equation of motion need be spec- 
ified, just the metric, to know how light is redshifted), 
this breaks down as soon as time dependence is explicit, 
e.g. by parametrizing q = qo + q\z. Thus no advantage 
exists for such a representation over the dynamical phase 
space. 



V. EXTENDING DYNAMICS 

We can now investigate whether the dynamics phase 
space w-w' is useful for physical theories beyond canon- 
ical, minimally coupled scalar fields. This includes for 
modified gravity or other theories where the quantities w 
and w' are effective quantities, defined in terms of the de- 
viation in the expansion rate from the matter dominated 
behavior, 



SH 2 = {H/H ) 2 ~n m a 
ld\n5H 2 
6 a ma 



(37) 
(38) 



possibly distinct from any physical pressure or dark en- 
ergy density. 

As already mentioned, phenomeno logical models such 
as 5H 2 ~ H a (jDvali fe Turnerl . 120031 ) fit within the freez- 
ing picture and the specific freezing region of the phase 
space, as illustrated in Fig. [3l Note that the case a = 1 
corresponds to th e dynamics of an extra dimensional 
braneworld model (iDeffavet. Dvali. fc Gabadadzei . l2002t 



iDvali. Gabadadze. fc Porrati L 2000[); suc h modelsare dis- 
cussed in more detail by Koyama ( 2007h in this volume. 

since the results of JEII] were discussed in terms of 
canonical, minimally coupled fields, let us examine the 
extension to noncanonical or coupled dark energy. 

k-essence: If we remove the canonical na- 
ture of the scalar field Lagrangian that involves 
an additive term linear in the kinetic energy, 
we have a class of theories known as k-e ssence 
(lArmendariz-Picon. Mukhanov. fc Steinhardtl |200(J ; 
IChiba. Okabe. fc Yamaguchil 12000) , with Lagrangians 
of the form 



C = V(4>)F(X), 



(39) 



where X = (d p 4>d tl (j))/2, i.e. in the absence of spatial in- 
homogeneities X is just the kinetic energy. Such models 
have some inspi rations from fi eld a nd str ing theory (for 
an overview see iNovello et all (e.e. 12005^ . can describe 
phantom fields with w < —1, can have sound speeds less 
than the speed of light (hence affecting structure forma- 
tion differently than quintessence) and can have attractor 
mechanisms to alleviate the fine tuning problem. 
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FIG. 3 Modifications to the Friedmann equation of the form 
H a lie in the freezing regime, despite possibly not arising 
from a simple scalar field. Moreover, they asymptotically ap- 
proach A along the lower boundary line w' = 3tu(l + w). The 
braneworld curve is shown solid to z = 0, with x's indicating 
z = 1,2,3. 



Without further specifying the functions V or F, it 
is difficult to say anything general about k-essence dy- 
namics. Purely kinetic k-essence, where V — constant, 
does have phase space trajectories limited to one side 
or the other of the line w' = 3w(l + w) correspond - 
ing to constant pressure (|Linderll200fibllSche rrerj,|2006). 
However kinetic k-essence can dynamically mimic (or be 
mimicked by) quintessence as long as the portion of the 
phase space trajectory of interest does not cross this line 
(Ide Putter fc Linden . [20071 ISenl 120061) . 

Coupled dark energy: The dark energy could in fact 
be not dark, that is it could interact non-gravitationally. 
From the dynamical perspective this creates an effective 
equation of state shifted from the bare one by the inter- 
action term, e.g. 



w cS = w 



r 

3H' 



(40) 



where T is the interaction appearing in the continuity 
equation 

p w = -3Hp w (l + w) + Tp Wl (41) 
representing a decay /creation process for ex ampl e. This 



representing a decay/ creation process tor ex ampl e, inis 
was set forth in early work by iTurnerl ( 19851 ) and iLinderl 
(1988a). Such coupling will shift the trajectories in the 



and t heir cosmological effects, have been considered; see , 
e.g.. (lAmendola fc Quercellinil . [200l iBarnes et all 120051: 
ILinderl l2005al) . However, concerns have been raised 
about the apparent strong effect of quantum corr ections 
on fields coupled to matter ( Doran fc Jackell [2002) . This 
can be avoided if one postulates that the potential con- 
sidered is really an effective low energy potential that 
just happens to take on a simple form as a result of com- 
plicated quantum loop corrections to the (in turn nec- 
ess arily complica ted) classical potential; see the article 
by iDurrerl (|2007l ) on low energy effective theories in this 
volume. 



Scalar-tensor gravity: Rather than coupling the 
dark energy to the matter sector of the La- 
grangian, one could make the coupling to grav- 
ity nonminimal. These are scalar-tensor theo- 
ries; see the article by iFrancaviglia fc Capoziellol 
(2007) in this volume. Coupling the quintessence 
field to the Ricci scalar, R/(8irG) -> F{(j))R 
i n the action, these extended quintess ence theories 
( Perrotta. Baccigalupi. fc Matarresell2000l ) can have var- 
ied dynamics depending on the form of F, along with 
an interesting attractor mechanism called the i?-boost 
([Baccigalupi. Matarrese. fc Perrottal 120001 ) . For a model 
with a cosmological constant potential, requiring consis- 
tency with solar system tests drives the equation of state 
very close to ic = —1 (within 10 -4 ) and with dynam- 
ics representative o f neit her freezing nor thawing fields 
dBaccigalupi et ail. 2007f). For a nother approach, see 
iNesseris fc Perivolaropoulosl ( 20071 ). 



Model Zoo: As the fertile imagination of children's au- 
thor Dr. Seuss envisioned an alphabet and animals "On 
Beyond Zebra" , so has the intense interest in the dark en- 
ergy mystery led to a zoo of models "On Beyond A" . The 
merest glimpse of a small fractio n of these includes: os- 
cillat i ng (see also slinky) models (iBarenboim fc Lvkkenl . 



w-w 1 phase space, allowing for dynamics outside the 
thawing and freezing regions. Many different couplings, 



120061 : IBarenboim. Mena Reaueio. fc Quiggj . [2006J) with 
dynamics corr esponding to a circle in phase space 
(ILinderl . l2006af) . mocker models that arc from matter 
like behavior to cosmological constan t like be havior along 
curves of w' — Cw(l + w) ([Linderl . l2006b[ ). closely re- 
lated to quartessence and Chaplygin gas models that 
attempt to unify dark mat t er an d dark energy (see 
iMakler. de Oliveira. fc Wagal (2003) for an overview), 
skating models that arc from free field behavior (w — 
+1), to cosmological constant like behavior along the 
curve w' = —3(1 — w 2 ), physically corr esponding to a 
field moving across a constant pot ential (jLinderl . [2005bJ; 
ISahlen. Liddle. fc Parkinson . 20051) (b u t also related to 
kinetic k-essence (Ide Putter fc Linderl 120071) ). and wet 
fluid (|Holman fc Naidul . 120041 ) (equivalent to the sum of 
a constant w compone nt and a cosm ological constant; 
cf. ClLD!) or leveling ([Linderl . |2006b) models that ap- 
proach a cosmological constant as the density nears a 
limiting value and have parabolic tracks - respectively 
w' = 3(1 + w)(w — W*) and w' — —3(1 + w)(w + w*). 
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VI. DYNAMICS AND GROWTH 



A. Thawing Physics 



The dynamics of the accelerating component affects 
the growth of structure in the universe through the ex- 
pansion rate. This provides a Hubble friction term op- 
posing gravitational instability (e.g. reducing the expo- 
nential Jeans growth in a static background to the power 
law growth in an expanding background). It also affects 
the matter source term Q m (a), i.e. the evolution of the 
homogeneous matter density, through the expansion, but 
to the extent that dark energy remains smooth on the rel- 
evant scales it does not directly source growth. Canon- 
ical scalar fields are very light, m < H, so t hey remain 
smoo th on scales less than the Hubble scale ( Ma et ail . 
1999). Therefore, within general relativity, the growth 
effects of dark energy follow directly from the expansion 
effects discussed in this article. A highly accurate fitting 
formula for the linear growth can be given in terms of 
fl m (a) and w(z = 1) through t he gravitational grow th 
index formalism (lLindeii[2005atlLinder fc Cahnl . I2007T) : 



9(a) 
7 



{5pm/ Pm)/a = 

0.55 + 0.05 [1 



e / o °(do/a)[n m (o)- 1, -l] / 42 x 

w(z = l)}, (43) 



is accurate to 0.2% compared to the numerical solution 
of the exact second order differential equation. Struc- 
ture formation in general requires treatment of fully non- 
linear growth through N-body numerical computations. 
Ea rly work with dyna m ical quintessence i nclude d that 
of iKlvpin etafl (|2003h : iLinder fc Jenkins! (|2003h . with 
many following investigations. 

When the physics of the cosmic acceleration has a grav- 
itational origin, or a dark energy component is not mini- 
mally coupled, additional terms enter into the growth, in- 
cluding new source terms such as from anisotropic stress 
and non-unity sound speed, and varying gravitational 
coupling. This breaking of the degeneracy between ex- 
pansion effects and growth effects offers a promising win- 
dow for identifying the fundamental physics, but is be- 
yond the sco pe of this article; see, e.g., the review by 
iLinderl (|2007ft for more details. 



Thawing models are defined by their departure from a 
cosmological constant-like state in the past to a dynam- 
ical, w =/= — 1, behavior today. This property of being 
frozen into a cosmological constant over much of the his- 
tory of the universe makes this class difficult to distin- 
guish from a cosmological constant without highly accu- 
rate cosmological data. Indeed, current observations are 
almost wholl y degenerate with the enti re thawing region 
as defined in ICaldwell fc Linden (l2005lh and if an effec- 
tive, constant w (e.g. a weighted average over the data 
sensitivity) is determined to equal —1 within 5% then 
we still have essentially no information on whether this 
is truly a cosmological constant A or any model in the 
entire half of the physical model space that is categorized 
as thawing. 

This challenge in uncovering the underlying physics 
makes this class useful as a testbed for the science reach 
of next generation experiments and for the role of phe- 
nomenological parameterization. We will particularly be 
interested in, of course, distinguishing thawing models 
from A and seeing dynamics such that w(z) ^ w cons t, 
but we also would learn physics more directly by verify- 
ing that the field started in a frozen state at early times 
and furthermore discerning its trajectory in phase space 
or at least its dynamical slop e parameter u>'/(l + w). 

Recall from ILinderl (|2006bl ) that 



= 3- 



2A + 3(l + w) 
1 - Y 



l + Y 



3w, 



(44) 
(45) 



where X = <fi/(H(j>) and Y = <p/V^,. Thus, constraining 
the dynamical slope parameter w' / (1 + w) directly leads 
to information about the field acceleration, friction, and 
potential tilt terms in the Klein-Gordon equation of mo- 
tion ©. 



B. Thawing Models 



VII. THAWING DARK ENERGY 

Let us now return to the fundamental mode picture 
of quintessential dynamics, presenting some new results 
on the specifics of determining the class of dark energy 
responsible for cosmic acceleration and the ability to zero 
in on characteristics within that class. 

While distinguishing the thawing class of dark energy 
from the freezing class would be a major accomplishment 
guiding us toward the fundamental physics behind dark 
energy, we can also examine thawing models in them- 
selves. These are among the best motivated physics, in- 
cluding radiatively stable PNGB pseudoscalar or axion 
models and familiar quadratic, quartic, and other renor- 
malizable potentials. 



We begin by examining three parameterizations of 
thawing fields, comparing their behavior and constraints. 
First is the standard parameterization of w(a) — wq + 
u> a (l — a), reviewed in £)IVI If we choose wo + w a = —1, 
then we see that at early times (a <C 1), this possessed 
w = —1. Furthermore, w' — —aw a so at early times 
w' = 0; thus this parameterization can describe a thaw- 
ing model. However, we have handcuffed this parameter- 
ization by doing this, reducing it to a single parameter, 
rather than a model with two degrees of freedom, putting 
it at a disadvantage. It is basically restricted to the tra- 
jectory w' = 1 + w. Nevertheless, we will see that it is 
able to describe reasonably most thawing models. The 
alternative is to retain the two parameters of wq , w a but 
at the price of not matching a cosmological constant at 
early times; since cosmological data weights the recent 
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universe more heavily, this is not a bad approximation. 
The energy density of w a models is 



Pw(a) 



a ~ 3(l + U)0+UJa)g — 3tU a (l— a) 



(46) 



The second parameterization w' = F(l + w) is moti- 
vated by PNGB mode ls, and is an excellent app roxima- 
tion to their dynamics ( Caldwell fe Linderl . l2005l ). with F 
inversely proportional to the symmetry b reaking energy 



scale /. For more on PN G B mo d els, see Fricma n et al 



1995); iKaloper fe Sorbol (|2006l ): iRosenfeld fe Frieman 



2006). These have fields starting frozen on their 1 



cos(0//) potential and, after the Hubble drag diminishes, 
are released to roll. Due to the change from convexity to 
concavity of the potential, they can have interesting dy- 
namics depending on the initial conditions. We assume 
they are not fine tuned in the sense of starting very near 
the top of their potential, nor have they already rolled 
through the minimum and ascended the potential. Even- 
tually the field will oscillate around the minimum (which 
looks quadratic, i.e. V ~ (j) n with n — 2, s o the effective 
equation of state w=(n- 2)/{n + 2) = (|Turnerl . [l983h 
as long as the oscillation period is short compared to the 
Hubble time), acting like dark matter before vanishing 
as the field comes to rest at zero potential. However, 
during the accelerating period, w' — F(l + w) accurately 
describes the dynamics; the equation of state has two 
parameters, the current equation of state wq and the dy- 
namical slope F, with 

1 + w = (1 + w ) a F . (47) 

The energy density of these models is 

Pw (a) = p w e^ 1+ ^ F ^- aF l (48) 

For the third parameterization, we craft a new model 
specifically following the physics of thawing, called the 
algebraic thawing model: 



1 + w = (1 + w ) aP 



1 



1 + ba 



l-p/3 

: , - (49) 



with parameters u>o, p (6 is fixed). Let us justify this 
form. As the field is released from the cosmological con- 
stant state, still in the matter dominated era t ~ a 3 / 2 , 
the dynami cs is given by X = 3/2, or w' = 3(l + w), as il- 
lustrated in ICaldwell fc Linderl (|2005l) . This implies that 
at early times 1 + w ~ a 3 . So far this is identical to the 
PNGB model with F — 3. To add some curvature into 
the trajectory in the phase plane w-w' , let us multiply 
this by a factor that bends the dynamics away from this 
line as the scalar field energy density becomes more im- 
portant, say Q w (a) q . In fact, to preserve the early time 
behavior, this factor must go to a constant at early times, 
so we use [a 3 n, w (a)) q , which is indeed constant at early 
times when w — ► — 1. The only problem with this is that 
the expression for the equation of state has become non- 
analytic. Even if we approximate il w (a) by some fixed 



function, say Q\ (a), then the equation of state is inter- 
twined with the present energy density parameter Q w , 
or f2 m , rather than being an independent quantity. For 
the final form we therefore replace the intruding density 
ratio - in this one place - with a constant b = 0.3. The 
equation of state is quite insensitive to this specific value, 
varying by less than 1% as b varies by 50%; of course the 
value of b is irrelevant as a — > 1 and for a<l. 
The dynamics of the algebraic thawing model is 



w' = (1 + w) 

and the energy density is 

"3(1 + w ) 



P 



1 + ba~ 3 



Pw (a) = p w exp 



ap 



[l-{aa 3 +0f /3 



(50) 



. (51) 



where a = 1/(1 + 6), = 6/(1 + 6). 

In a clever analysis, ICrittenden. Maierotto. fc Piazzal 
(120071) came up with a similar model by analyzing a slow 
roll-like field expansion, assuming a particular combina- 
tion V,f,/[V(1 + X/3)] can be Taylor expanded about the 
present. After some approximations they take 1 + w ~ 
a p d,\(a) 1 ~ p ^ 3 . However, this form still entangles w and 
the present matter density, and in fact a more exact so- 
lution of the field expansion equations works worse! The 
basic problem is that even for thawing fields there is no 
reasonable slow roll or field expansion approximation. 
Even for their less extreme model with wo = —0.8, the 
field still traverses A0 ~ 0.4Mp. The algebraic form 
Eq. ([49]) in fact gives more accurate equations of state 
for the cases they illustrate. 



C. Discriminating Thawing 

We can now use the w a , PNGB, and algebraic models 
to examine the constraints, and parameter dependence 
of the constraints, from future data on the dynamics of 
quintessence. For each model we have two equation of 
state fit parameters: (wo,w a ), (wo,F), or (wo,s), where 
s = w' /(l+Wo) is the dynamical slope at present (just as 
F is the dynamical slope, constant for all times). From 
the estimation of these parameters (marginalizing over 
the matter density, in a flat universe, and other param- 
eters such as the supernovae absolute magnitude), and 
their covariances, we can find the constraints on w and 
w' at any redshift, giving confidence contours in the w-w' 
phase plane. 

For future data we consider Type la supernovae dis- 
tances from z = — 1 . 7, wit h systematics, of SNAP qual- 
ity (see, e.g., SNAP (2004)), plus the reduced distance 
to the CMB last scattering surface, of Planck quality 
(0.7% fractional precision). The fiducial cosmology has 
wo = —0.9 and present dynamical slope 1.5, and likeli- 
hoods are approximated as Gaussians in a Fisher infor- 
mation analysis. The CMB data in fact has little lever- 
age on the equation of state, because for all the thawing 
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models the high redshift equation of state goes to a cos- 
mological constant. We have checked that adding baryon 
acoustic oscillation angular distance measurements at 1% 
precision or a matter density prior of 0.005 (roughly mim- 
icking weak gravitational lcnsing constraints) does little 
to improve the constraints. 

Figures show the w-w' constraints for the three 
models at four redshifts. We exhibit the 68% confidence 
level contours at z = 0, at the redshift where w and w' 
are decorrelated, giving vertical/horizontal ellipses, and 
at high redshift, z > 1. The phase space trajectory is 
marked by the x's at each of the four redshifts. Note that 
the confidence contours vary between the models, espe- 
cially when evaluated at the present, and this may lead to 
concerns about parameterization dependence. However, 
as we will see, the qualitative answers to the important 
physical questions remain independent of the parameter- 
ization. 
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FIG. 4 Together with Figs. EUl this fi gure for the w a thawing 
model illustrates constraints on the dynamical behavior of 
three thawing models at four redshift snapshots. While the 
z = behavior is poorly limited by future data, taking into 
account the dynamical history still allows distinction of the 
fiducial wo — —0.9, w' = 0.15 model from a cosmological 
constant and from the freezing class of physics. 



While constraints on the present dynamical state, i.e. 
Wo and w' , are relatively weak, in each of the parameter- 
izations they are still sufficient to distinguish the fiducial 
model w — —0.9, w' — 0.15 from a cosmological con- 
stant. (Note this is despite the uncertainty a(wo) ~ 0.14 
from the algebraic thawer, the weakest model at z = 
- one must take into account the contour orientation in 
the phase plane.) At z = 0, however, the models can- 
not distinguish thawing from freezing, or from a constant 
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As Fig. H for the PNGB thawing model. 
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FIG. 6 As Fig. [4] for the algebraic thawing model. 



equation of state w C onst = —0.9. Using the information 
from throughout the dynamical history greatly improves 
the situation. At some redshift, z « 0.2 — 0.3 in the 
cases here, the dynamical variables w and w' decorre- 
late and the contours become vertical. This gives the 
greatest distance between the constraint contour and the 
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cosmological constant, showing clear distinction, and the 
intersection of the ellipse with the w' = axis also pro- 
vides the minimal variance estimate on the instantaneous 
equation of state value. Such a decorrelation redshift is 
sometimes called a pivot redshift. Gcnerically there can 
be more than one decorrelation redshift, and for the mod- 
els where w is not a linear function of the parameters 
we exhibit the contours at the second of these redshifts, 
z w 0.7 — 1.4. This provides a minimum variance esti- 
mate of the instantaneous time variation of the equation 
of state. 

Note that the confidence contours at each redshift are 
distinct from the cosmological constant, showing that 
future data can distinguish thawing models from A (at 
least at la for this fiducial cosmology). Furthermore, the 
early time contours (except in the w a case) distinguish 
the thawing model from models with constant equation 
of state, thus exhibiting the presence of dynamics. The 
early time contours also draw away from the freezing re- 
gion of the phase plane, so the data can indeed guide us 
to the correct class of physical origin. These are all im- 
portant physical insights that are not parameterization 
dependent. Gains are more modest in zeroing in on a 
specific thawing model and these are more sensitive to 
parameterization. At early times, the form of the al- 
gebraic thawer forces the contour to prefer a dynamical 
slope near 3. However the PNGB and w a cases do not 
impose such preferences since the slope is a free fit pa- 
rameter. They do constrain w' / (1 + w) to a subset of the 
thawing region, rather than the full range of 1-3. 

It is heartening that the physical insights can be ex- 
pected to be as clear as indicated, and not particularly 
dependent on the specific parameterization. The issue of 
fitting the dynamical behavior of dark energy (especially 
when restricted to two parameters, as seems likely from 
realistic next generation data accuracy), is a fascinating 
one. Use of a global parameterization like (wo,w a ) al- 
lows a good fit for models over the whole phase plane, 
but one can imagine that as we close in on the physical 
origin of dark energy, e.g. narrowing in on thawing mod- 
els, we may move to more specific parameterizations such 
as the algebraic thawing model. On the other hand, per- 
haps specific physical benchmark models, such as PNGB 
or motivated scalar field potentials, will then be of most 
use. 



VIII. CONCLUSION 

Dynamics, of quintessence and of the accelerating 
physics in general, can provide considerable insight into 
the nature of the new component or new physical law 
dominating our present universe. Fundamental modes of 
the physics lead to well defined, distinct regions of w-w' 
phase space that next generation cosmological probes will 
be able to test and distinguish. Just as we build our phys- 
ical intuition in early universe inflation with single field 
models leading to consistency relations, the fundamental 



modes of dark energy - the quintessence of dynamics - 
are a useful foundation. 

Model independent parametrization, with a strong 
physical basis, plays an important role, even if stringent 
constraints will be limited to two parameters such as the 
tilt from a cosmological constant, 1 + w, and a varia- 
tion w'. Nevertheless, this is as much as we expect from 
inflation as well, while for dark energy we have added 
complications due to the incomplete dominance of dark 
energy. 

Sensitivity to dynamics is a requirement to make 
progress in understanding the nature of cosmic acceler- 
ation. Once we begin to zero in on a class of physics, 
model independence may give way to specific discrimi- 
nating approaches such as the thawing analysis presented 
here. Models for the equation of state which depend non- 
linearly on the time variation parameter also possess min- 
imum variance, or pivot, redshifts for the time variation, 
z p i , and this may prove a useful tool. 

Dynamics alone, whether by its characterization or ab- 
sence, will not fully solve the dark energy enigma. The 
cosmic expansion history must be properly compared 
with the cosmic growth history to reveal extensions to 
gravitational physics or microphysics. We have scarcely 
addressed this important subject here, nor have we said 
why in the presence of dynamics A should not still exist, 
at a much larger energy density than the present, causing 
an abnegation of the universe we observe. 

Ten years passed from the time the basic physics and 
cosmology for the accelerating universe were in place un- 
til the first convincing observational evidence for its re- 
ality; since then another ten years of work on all fronts 
have passed. There is clearly still an enormous amount 
of exciting and challenging work ahead, and the an- 
swers, whatever they are and whenever they come, will 
revolutionize our understanding of gravitation, quantum 
physics, cosmology, and the fate of our universe. 
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